R Notebook

This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Reverse Engineer 2 Notebook

Objective: Analyze data from the Blair and Elrich campaigns for county executive from 2018 and 2022.

###Question 1:

First thing we did is load up the sets and clean the data:

##Blair 2018 contribs:

blair_contribs_2018 <-read.csv("data2/blair_contribs_2018.csv")

blair_contribs_2018_cleaned <-blair_contribs_2018 %>% 
clean_names  
  
  
blair_contribs_2018_cleaned %>%   
mutate(contribution_date = mdy(contribution_date)) 

blair_contribs_2018_cleaned

###Blair 2022 contribs

blair_contribs_2022 <-read.csv("data2/blair_contribs_2022.csv")

blair_contribs_2022_cleaned <- blair_contribs_2022 %>% 
clean_names 

blair_contribs_2022_cleaned %>% 
mutate(contribution_date = mdy(contribution_date)) 

blair_contribs_2022_cleaned

###Blair 2022 expenses

blair_expenses_2022 <-read.csv("data2/blair_expenses_2022.csv")

blair_expenses_2022_cleaned <-blair_expenses_2022 %>% 
  clean_names 

blair_expenses_2022_cleaned %>% 
mutate(contribution_date = mdy(expenditure_date)) 

blair_expenses_2022_cleaned

###Blair 2018 expenses

blair_expenses_2018 <-read.csv("data2/blair_expenses_2018.csv")

blair_expenses_2018_cleaned <-blair_expenses_2018 %>% 
  clean_names()

blair_expenses_2018_cleaned %>% 
  mutate(contribution_date = mdy(expenditure_date)) 

blair_expenses_2018_cleaned 

##Elrich 2018 contribs:

elrich_contribs_2018 <-read.csv("data2/elrich_contribs_2018.csv")

elrich_contribs_2018_cleaned <-elrich_contribs_2018 %>% 
  clean_names()

elrich_contribs_2018_cleaned %>% 
  mutate(contribution_date = mdy(contribution_date)) 

elrich_contribs_2018_cleaned

###Elrich Contribs 2022

elrich_contribs_2022 <-read.csv("data2/elrich_contribs_2022.csv")

elrich_contribs_2022_cleaned <-elrich_contribs_2022 %>%
  clean_names()

elrich_contribs_2022_cleaned %>% 
  mutate(contribution_date = mdy(contribution_date)) 

elrich_contribs_2022_cleaned

###Elrich 2018

elrich_expenses_2018 <-read.csv("data2/elrich_expenses_2018.csv")

elrich_expenses_2018_cleaned <-elrich_expenses_2018 %>% 
  clean_names()

elrich_expenses_2018_cleaned %>% 
    mutate(contribution_date = mdy(expenditure_date)) 

elrich_expenses_2018_cleaned

###Elrich Expenses 2022

elrich_expenses_2022 <-read.csv("data2/elrich_expenses_2022.csv")

elrich_expenses_2022_cleaned <-elrich_expenses_2022 %>% 
  clean_names()

elrich_expenses_2022_cleaned %>% 
  mutate(contribution_date = mdy(expenditure_date)) 

elrich_expenses_2022_cleaned

##County primary election results:

dem_precincts_18 <- read.csv("data2/dem_precincts_2018.csv")

dem_precincts_22 <- read.csv("data2/dem_precincts_2022.csv")

dem_county_22 <- read.csv("data2/dem_county_2022.csv")

dem_county_18 <- read.csv("data2/dem_county_2018.csv")

###Question 1 1. How much money did David Blair and Marc Elrich fund themselves in 2022 vs. 2018? Blair is a businessman and millionaire, so it makes sense that he’s funding his own campaign. But how much exactly is he putting into it, and how does it compare to the past election and to how much Elrich’s campaign is raising?

###2022:

blair_contribs_2022_cleaned %>%
filter(contributor_name == "BLAIR  DAVID  THOMAS") %>%
group_by(contributor_name) %>%
summarize(total_blair = sum(contribution_amount))

elrich_contribs_2022_cleaned %>%
group_by(contributor_name) %>%
summarize(total_elrich = sum(contribution_amount)) %>%
arrange(desc(total_elrich))

###2018:

blair_contribs_2018_cleaned %>%
filter(contributor_name == "BLAIR  DAVID  THOMAS") %>%
group_by(contributor_name) %>%
summarize(total_blair = sum(contribution_amount))

elrich_contribs_2018_cleaned %>%
group_by(contributor_name) %>%
summarize(total_elrich = sum(contribution_amount)) %>%
arrange(desc(total_elrich))

###Question 2: Who were the top 5 contributors to Blair and Elrich in 2022? How about in 2018? What are their connections to the candidates?

###2022

top_blair_contribs_22 <- blair_contribs_2018_cleaned %>%
  mutate(contributor_name = case_when(
    contributor_name == "BLAIR  DAVID  THOMAS" ~ "Blair  David  Thomas",
    TRUE ~ contributor_name
  )) %>%
  group_by(contributor_name) %>%
  summarize(total = sum(contribution_amount)) %>%
  arrange(desc(total)) %>%
  head(5)
top_blair_contribs_22

top_elrich_contribs_22 <- elrich_contribs_2022_cleaned %>%
  group_by(contributor_name) %>%
  summarize(total = sum(contribution_amount)) %>%
  arrange(desc(total)) %>%
  head(6)
top_elrich_contribs_22

###2018

top_blair_contribs_18 <- blair_contribs_2018_cleaned %>%
  mutate(contributor_name = case_when(
    contributor_name == "BLAIR  DAVID  THOMAS" ~ "Blair  David  Thomas",
    TRUE ~ contributor_name
  )) %>%
  group_by(contributor_name) %>%
  summarize(total = sum(contribution_amount)) %>%
  arrange(desc(total)) %>%
  head(5)
top_blair_contribs_18

top_elrich_contribs_18 <- elrich_contribs_2018_cleaned %>%
  group_by(contributor_name) %>%
  summarize(total = sum(contribution_amount)) %>%
  arrange(desc(total)) %>%
  head(6)
top_elrich_contribs_18

###Question 3: David Blair got more early/election day votes in 2022. Elrich got more Mail-in votes in 2022. How did that compare to 2018? People want to know how Marc Elrich won both times, first by 72 voters in 2018 and then 32 in 2022. Where did Marc do well in the three categories: Early voting, election day, and mail-in votes. Did this sway the result at all?

dem_county_18_cleaned <-dem_county_18 %>% 
  clean_names()

dem_county_22_cleaned <-dem_county_22 %>% 
  clean_names()

blair <- dem_county_18_cleaned %>%
filter(candidate_name == "David Blair")
elrich <-  dem_county_18_cleaned %>%
  filter(candidate_name == "Marc Elrich")
blair_elrich <-  bind_rows(blair, elrich)
blair_elrich

###Question 4: Which parts of the county voted for Elrich and which voted for Blair based on precinct-level voting? Are there differences in demographics of those areas? What about income?

\

library(sf)

Linking to GEOS 3.9.1, GDAL 3.4.3, PROJ 7.2.1; sf_use_s2() is TRUE

moco_prec_2022 <- st_read("data2/MontMD_2022/MontMD_2022.shp") %>% 
  st_zm(drop=TRUE)

Reading layer `MontMD_2022' from data source 
  `C:\Users\merca\Documents\GitHub\data_journalism_fall_2022\major_assignments\reverse_engineering\data2\MontMD_2022\MontMD_2022.shp' 
  using driver `ESRI Shapefile'
Simple feature collection with 261 features and 5 fields
Geometry type: MULTIPOLYGON
Dimension:     XYZM
Bounding box:  xmin: 1162792 ymin: 461729.6 xmax: 1344207 ymax: 614681.3
z_range:       zmin: 0 zmax: 0
m_range:       mmin: 0 mmax: 0
Projected CRS: NAD83 / Maryland (ftUS)

glimpse(moco_prec_2022)

Rows: 261
Columns: 6
$ NAME     <chr> "MONTGOMERY PRECINCT 07-032", "MONTGOMERY PRECINCT 09-015", "MONTGOMERY PRECINCT 13-056", "M…
$ NUMBER   <chr> "07-032", "09-015", "13-056", "13-016", "01-006", "04-027", "08-008", "08-016", "10-001", "0…
$ JURSCODE <chr> "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MONT", "MON…
$ VOTESPRE <chr> "007-032", "009-015", "013-056", "013-016", "001-006", "004-027", "008-008", "008-016", "010…
$ COUNCIL  <chr> "1", "3", "6", "4", "7", "6", "7", "7", "1", "6", "7", "1", "1", "7", "4", "3", "3", "2", "6…
$ geometry <MULTIPOLYGON [US_survey_foot]> MULTIPOLYGON (((1289150 485..., MULTIPOLYGON (((1251625 546..., MU…

moco_boundaries <- st_read("data2/moco_boundary.gdb")

Reading layer `CNTY_BNDY' from data source 
  `C:\Users\merca\Documents\GitHub\data_journalism_fall_2022\major_assignments\reverse_engineering\data2\moco_boundary.gdb' 
  using driver `OpenFileGDB'
Simple feature collection with 1 feature and 3 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: -77.52769 ymin: 38.93425 xmax: -76.88764 ymax: 39.35434
Geodetic CRS:  WGS 84

glimpse(moco_boundaries)

Rows: 1
Columns: 4
$ NAME         <chr> "MONTGOMERY COUNTY"
$ SHAPE_Length <dbl> 1.895903
$ SHAPE_Area   <dbl> 0.1367927
$ SHAPE        <MULTIPOLYGON [°]> MULTIPOLYGON (((-77.18523 3...

moco_boundaries %>%
  ggplot() +
  geom_sf() +
  theme_minimal()

moco_prec_2022

Simple feature collection with 261 features and 5 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 1162792 ymin: 461729.6 xmax: 1344207 ymax: 614681.3
Projected CRS: NAD83 / Maryland (ftUS)
First 10 features:
                         NAME NUMBER JURSCODE VOTESPRE COUNCIL                       geometry
1  MONTGOMERY PRECINCT 07-032 07-032     MONT  007-032       1 MULTIPOLYGON (((1289150 485...
2  MONTGOMERY PRECINCT 09-015 09-015     MONT  009-015       3 MULTIPOLYGON (((1251625 546...
3  MONTGOMERY PRECINCT 13-056 13-056     MONT  013-056       6 MULTIPOLYGON (((1301382 519...
4  MONTGOMERY PRECINCT 13-016 13-016     MONT  013-016       4 MULTIPOLYGON (((1298182 487...
5  MONTGOMERY PRECINCT 01-006 01-006     MONT  001-006       7 MULTIPOLYGON (((1268938 556...
6  MONTGOMERY PRECINCT 04-027 04-027     MONT  004-027       6 MULTIPOLYGON (((1280700 512...
7  MONTGOMERY PRECINCT 08-008 08-008     MONT  008-008       7 MULTIPOLYGON (((1282944 540...
8  MONTGOMERY PRECINCT 08-016 08-016     MONT  008-016       7 MULTIPOLYGON (((1295564 548...
9  MONTGOMERY PRECINCT 10-001 10-001     MONT  010-001       1 MULTIPOLYGON (((1247972 509...
10 MONTGOMERY PRECINCT 04-026 04-026     MONT  004-026       6 MULTIPOLYGON (((1285712 505...

head(moco_prec_2022)

Simple feature collection with 6 features and 5 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 1246398 ymin: 481181.4 xmax: 1302601 ymax: 556657.3
Projected CRS: NAD83 / Maryland (ftUS)
                        NAME NUMBER JURSCODE VOTESPRE COUNCIL                       geometry
1 MONTGOMERY PRECINCT 07-032 07-032     MONT  007-032       1 MULTIPOLYGON (((1289150 485...
2 MONTGOMERY PRECINCT 09-015 09-015     MONT  009-015       3 MULTIPOLYGON (((1251625 546...
3 MONTGOMERY PRECINCT 13-056 13-056     MONT  013-056       6 MULTIPOLYGON (((1301382 519...
4 MONTGOMERY PRECINCT 13-016 13-016     MONT  013-016       4 MULTIPOLYGON (((1298182 487...
5 MONTGOMERY PRECINCT 01-006 01-006     MONT  001-006       7 MULTIPOLYGON (((1268938 556...
6 MONTGOMERY PRECINCT 04-027 04-027     MONT  004-027       6 MULTIPOLYGON (((1280700 512...

write_csv(moco_prec_2022, "data2/moco_prec_2022.csv")

moco_prec_2022 %>%
  ggplot() +
  geom_sf() +
  theme_minimal()

moco_22_pre_results<-read.csv("data2/Moco_22_pre_results.csv")
head(moco_22_pre_results)

moco_22_pre_results_cleaned <-moco_22_pre_results %>% 
  clean_names()
moco_22_pre_results_cleaned

moco_22_pre_results_filtered<-moco_22_pre_results_cleaned %>%
filter(candidate_name == "David T. Blair" | candidate_name == "Marc Elrich")
moco_22_pre_results_filtered

blair_results_2022_filter<-moco_22_pre_results_filtered %>% 
  filter(candidate_name == "David T. Blair")
blair_results_2022_filter

blair_results_2022_filter_total <-blair_results_2022_filter%>% 
  mutate(total_votes_blair = early_votes + election_night_votes + mail_in_ballot_1_votes + provisional_votes + mail_in_ballot_2_votes)
blair_results_2022_filter_total

blair_results_joined_2022<- inner_join(moco_prec_2022, blair_results_2022_filter_total, by=c("VOTESPRE"="election_district_precinct"))

###Blair Map####

ggplot() +
  geom_sf(data=blair_results_joined_2022, aes(fill=total_votes_blair)) +
  scale_fill_viridis_b(option="magma")+
  theme_minimal()

####Elirich map####

moco_22_pre_results_filtered_elrich <-moco_22_pre_results_cleaned %>%
filter(candidate_name == "Marc Elrich")
moco_22_pre_results_filtered_elrich

elrich_22_filtered_total <-moco_22_pre_results_filtered_elrich%>% 
  mutate(total_votes_elrich = early_votes + election_night_votes + mail_in_ballot_1_votes + provisional_votes + mail_in_ballot_2_votes)
elrich_22_filtered_total

elrich_results_joined_2022_map<- inner_join(moco_prec_2022, elrich_22_filtered_total, by=c("VOTESPRE"="election_district_precinct"))

ggplot() +
  geom_sf(data=elrich_results_joined_2022_map, aes(fill=total_votes_elrich)) +
  scale_fill_viridis_b(option="magma")+
  theme_minimal()

###percentage_change***

blair_results_2022_filter_total

elrich_22_filtered_total

both_joined_2022<- left_join(blair_results_2022_filter_total, elrich_22_filtered_total, by=c("county"="county"))

both_joined_2022

both_joined_mutated<-both_joined_2022 %>% 
  mutate(difference = total_votes_elrich-total_votes_blair)
both_joined_mutated

NA

both_joined_map<- inner_join(moco_prec_2022, both_joined_mutated, by=c("VOTESPRE"="election_district_precinct.x"))

ggplot() +
  geom_sf(data=both_joined_map, aes(fill=difference)) +
  scale_fill_viridis_b(option="magma")+
  theme_minimal()

library(ggplot2)

ggplot(results_joined_2022)+
  geom_map(aes(map_id = region))+
  expand_limits()

###Question 5: Where did David Blair and Marc Elrich spend their money on campaign finance resources? Ie: Meta advertisements, yard signs, TV, consulting? This would be interesting to know since Blair made a very big concerted effort to build his profile during the campaign while Elrich relied on incumbency and mainly TV ads to help him.

blair_expenses_2022_cleaned %>%
group_by(expense_category) %>%
summarize(total_category = sum(amount)) %>%
arrange(desc(total_category))
#spent nearly $3.1 mil in 2022 ... top 3 are Media, Salarties and Direct Mail by Mail House (R)

#dive into media

blair_expenses_2022_cleaned %>%
filter(expense_category == "Media") %>%
group_by(expense_purpose) %>%
summarize(total_category = sum(amount)) %>%
arrange(desc(total_category))

Spent 1773.00 on online advertising

elrich_expenses_2022_cleaned %>%
group_by(expense_category) %>%
summarize(total_category = sum(amount)) %>%
arrange(desc(total_category))

Elrich spent $528393.29 on media which is significantly less than Blair

elrich_expenses_2022_cleaned %>%
filter(expense_category == "Media") %>%
group_by(expense_purpose) %>%
summarize(total_category = sum(amount)) %>%
arrange(desc(total_category))

#Elrich spent 32K on consulting fees.

This is what I wrote for No. 5 “There are signifant differences in spending between Elrich and Blair here and it shows how the capital on hand can really influence the means a candidate has to win or lose a race. Elrich in 2022 went all-in on TV spending. Blair spent a lot more in 2018 as well.”